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ABSTRACT 

We present a new probe of primordial non-Gaussianity via the 21cm radiation from 
minihalos at high redshifts. We calculate the fluctuations in the brightness temperature 
(measured against the cosmic microwave background) of the 21cm background from 
minihalos containing HI at redshift ~ 6 — 20, and find a significant enhancement due to 
small non-Gaussianity with amplitude / NL < 1. This enhancement can be attributed 
to the nonlinear bias which is strongly increased in the presence of non-Gaussianity. 
We show that our results are robust against changes in the assumed mass function 
and some physical aspects of minihalo formation, but are nevertheless sensitive to 
the presence of strong radiation sources within or around the minihalos. Our findings 
are relevant for constraining and searching for small primordial non-Gaussianity with 
upcoming radio telescopes such as LOFAR and SKA. 



1 INTRODUCTION 



In this era of "precision cosmology" , many open questions in 
cosmology will be addressed with the highly anticipated re- 
sults from the Planck satellit^ and a host of other ambitious 
ground and space-based experiments. One of the key goals of 
these experiments is to establish the statistical nature of the 
primordial cosmological perturbations, which, according to 
the single-field, slow-roll model o f inflation, should be very 

Maldacenal <|2003h . for re- 



close to a Gaussian field (see e .g. ^_ 

views see lBartolo et all (|2004h : IChenl (|2010h ). A detection of 



any primordial non-Gaussianity (NG) would therefore sig- 
nify a deviation from the simplest model of inflation and 
hint at new physics in the early Universe. 

At present, the most stringent constraint on NG comes 
from measurements of the anisotropies in the cosmic mi- 
crowave background (CMB). In the "local" model of NG 
parametrized by a constant amplitude, / nl, the WMAP 
satell ite reported a limit /nl = 32 ±42 (2a) (Ko matsu et al.l 
I'JUI il l whilst Planck is expected to improve this limit to 
I/nl| < 5. 

On smaller physical scales, galaxy clusters have been 
shown to be effective in constraining NG via the number 
counts of rare obje c ts and measurement s of th e bias (see e.g. 
IScoccimarro et al.1 (12004); I Dalai et alJ il200sl); ISlosar et al. 



Desiacques fc Seliakl (|201Ch : 
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A positive /nl, for instance, will increase the num- 
ber of galaxy clusters at high redshifts and induce scale- 
dependence in the bias. Realistically, however, uncertain- 
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ties in the mass-observable relation, redshift determination 
and other systematics will most likely li mit the constrainin g 
power of cluster surveys to |/nl| > 10. ISlosar" et all (|200Sh . 
for example, obtained the 2a limit of —29 < /nl < 70 from 
the analysis of various large-scale-structure datasets. 

More recently, it has been shown that NG leaves im- 
prints in the 21cm radiation due to spin-flip transitions 
of neutral hydrogen during the epoch of reionization. This 
transition occurs at an excitation temperature T* = 68 mK 
and rest-frame frequency vq — 1.42 GHz (now redshifted 
to the radio band). These 21cm signals are expected to be 
measured over a wide range of redshifts by radio facilities 
such as LOFA ffl MWA0 and SKA0. It was first shown in 
ICooravl (|2006h that an ideal measurement of the 3-point cor- 
relation in the 21cm fluctuations from z ~ 50 can, in prin- 
ciple, pr obe /nl as small as 10 -2 (see also iPillepich et al.l 
l|2007h h IJoudaki et all (|201lh considered a more conserva- 
tive approach of measuring the 21cm power spectrum and 
concluded that |/nl| ~ 1 is within reach of th e next gener- 
ation of radio telescopes. iTashiro fc Sugivamal ((2012) anal- 
ysed the number counts of ionized bubbles in 21cm maps and 
found them to be se nsitive to |/nl| ~ 100. Most recently, 
ITashiro fc Hoi (|2012o argued that the correlation between 
the CMB and 21cm fluctuations will be capable of detecting 
|/nl| ~ 100 and will be useful for removing foregrounds and 
systematics. 

In this work, we present a new connection between NG 



2 Low- Frequency Array (www.lofar.org) 

3 Murchison Widefield Array (www.mwatelescope.org) 

4 Square-Kilometre Array (www.skatelescope.org) 



2 Chongchitnan and Silk 



and 21cm radiation, namely, via the 21cm signals from mini- 
halos (MHs), which are some of the earliest cosmic struc- 
tures to have formed post-recombination. MHs are typi- 
cally virialised clouds with mass 10 4 - 1O 8 M , compris- 
ing dark and baryonic matter, which have cooled to be- 
low ~ 10 4 K. These MHs host such a high density of neu- 
tral hydrogen that their 21cm signal can appear in emis- 
sion or absorption against the CMB, result i ng in a 21cm 
"forest" of spectral lines (llliev et al.l 20021; IShapiro et al.l 
l2006tlFurlanetto fc Loebll2002l ; lMeiksinll201ll ). Through the 
changes in number density and bias, we will show that the 
fluctuations in the 21cm emission from MHs provides a new 
probe of small NG with /nl J$ 1 and is potentially detectable 
by LOFAR and SKA. 

Throughout this work, we assume a flat ACDM cos- 
mology with matter density Q m = 0.276, baryon den- 
sity Q b = 0.046, Hubble constant h — 0.7, spectral index 
n s — 0.96 and power-spectrum normalisation erg = 0.81. 



2 21CM EMISSION FROM MINIHALOS 



We fo l low the basic treatme nt of halos in IShapiro et al.l 
l|l999l ): llliev fc Shapiro! l|200ll ) where a MH of a given mass 
is modelled as a "truncated isothermal sphere" with radius 
rt, temperature Tk, dark- matter velocity dispersion av, and 
density profile p(r). 

We begin by considering a single MH. Its 21cm signal 
may appear in emission or absorption against the CMB de- 
pending on the spin temperature, 7s, which is determined by 
1) energy exchanges between Hi-bound electrons and CMB 
photons, 2) collisions between atoms, 3) interactions be- 
tween electrons and Lym an-a (Lyq) photon s through th e 
Wouthuysen-Field effect <|Wouthuvsenl 1 19521 ; iFieldl 1 19581 ). 
We can write T s = (T C mb + y a T a + y c T K )/(l + y a + y c ), 
where y a , y c are the Lyq and collisional coup ling constants 
|Madaii et al.ll 19971 ; lAUison fc Dalgarnoll 19691 ) and T a ^ T K 
is the Lya colour temperature. Assuming for now that bright 
UV and X-ray sources have yet to form, or that the MH is 
isolated from such sources, we can set y a = 0. 

The observed brightness temperature at comoving dis- 
tance r from the centre of the halo is given by 



T b (r) = T C MB(0)e 



-r(r) 



+ / T s e- T dr\ 



(1) 



where Tcmb(^) = 2.73(1 + z) K. The optical depth, r, 
of neutral hydrogen to photons at rest-frame frequency, v, 
along a line of sight with impact parameter, a (in unit 
of rt,) from the centre of the MH can be expressed as 
jFurlanetto fc Loeblliool ) 



T(U) 



3c 2 A 10 T t 
32ttj/2 



n m {l)<t>{v,l) 



dR, 



(2) 



where A 10 = 2.85 x 10" 15 s" 1 and 7? and I are radial comov- 
ing distances satisfying £ 2 = R 2 + (art) 2 , with R — at the 
centre of the MH. The number density of neutral hydrogen 
in the MH is nm ~ (1 — Y)(Q b /Q m )(p/mH), where Y is the 
helium fraction and ran is the mass of a hydrogen atom. 
The intrinsic line profile <j>(v) is modelled as a Doppler- 
broadened form <j)(v) = (Au^/n) -1 exp (— \{y — i/q)/Av\ 2 ), 
with Ais — (uo I c) ^/2kBTK (-?) / tuh , and fcs the Boltzmann 
constant. 



In the special case when the line is unbroadened, 4>{v) = 
5(y — uq), the optical depth red uces to that of the unper- 
turbed IGM patch at redshift z (|Madau et alJll997h 



tigmU) 



3c 3 A 10 T t n m (z) 
32tt^T s (z)H(z) 



(3) 



Using (Hll-([3ll, we can rewrite the brightness temperature as 

Ts{£)e -T(»,R)dL dRj (4) 

- R nm(e')<Kv,e') 



T b {v) = TcM B e- TM + 



t(u,R) = tigm + 



3c 2 A w n 
32tvu 2 



Ts(t>) 



-dR'. (5) 



Finally, the observed 21cm brightness temperature of a sin- 
gle MH with respect to the CMB is given by 



ST b = -^--Tcmb(O), 

1+2 



(6) 



where T b is averaged over the halo cross-section A = irr 2 . 

If we now consider an ensemble of MHs in the mass 
range [M m i n , M max l, the mean 21cm emission from an en- 
semble is given by <|lliev et alj|2003 ) 



ST h = 



c(l + z) 4 
v H(z) 



fill 

Au cB ST b (M)A m dM. (7) 



where Av c g = [cj>(vo)(l + z)]^ 1 is the effective redshifted 
linewidth. Various prescriptions for the mass function, 
dn/dM, will be compared later on. We set M max to cor- 
respond to the virial temperature of 10 4 K whilst M m i n is 
set by the Jeans mass, Mj. 

The key observable relevant for the upcoming radio ar- 
rays is the rms fluctuations in the 21cm emission. For a 
pencil-beam survey with frequency bandwidth and an- 
gular size A9, the amplitude of the 3a fluctuation is 

(5T 6 2 ) 1/2 = 3a p (Au, A6)f3(z)5T bl (8) 

(see e.g. iDodelsonl (2003) for the calculation of the variance 
<T P in a cylinder). Here, P(z) is the weighted average of the 
bias b(M,z), defined as the ratio of the 2-point correlation 
for density peaks corresponding to a MH of mass M, and 
that of dark-matter density fluctuations (detail in the next 
section). 



5™—b{M,z)Hm) 

Jk£ mm 



T dM 



An in 



Hm)^jdM 



(9) 



/m. 

where .F(m) oc T b r 2 av is the effective flux from the MHs. 
3 EFFECTS OF /nl < 1 

If Planck rules out |/nl| > a few, then it would appear 
extremely difficult for large-scale-structure probes to ever 
improve on, or even corroborate, /nl constr aints from the 
CMB (unless /nl is fc-dependent, see e.g. lLoVerde et al.l 
(2008)). Unlike galaxy clusters, MHs are neither very rare 
nor very massive, so naively we expect the enhancement in 
their number counts from /nl < 1 to be undetectably small. 
However, at redshift > 6, these MHs are strongly biased non- 
linear objects and the scale-dependent effects on the bias can 
be much more dramatic t han those on clus ters. 

In the Gaussian case, llliev et all (120031) used the nonlin- 
ear bias approach of IScannapieco fc Barkanal (|2002T > (based 



21cm radiation from minihalos as a probe of small non-Gaussianity 3 



on excursion set theory) to study 21cm emission from MHs. 
They found very good agreement for (<5T 2 ) between analytic 
prediction and simulation and concluded that, for all prac- 
ti cal purposes, one cou ld, for instance, use the bias obtained 
in I Mo fe White! (1 19961 ) (using the peak-background split ap- 
proach). In this work, we extend this line of investigation to 
non-Gaussian scenarios. 

A number of previous authors have investigated the 
effect of non-Gaussianity on the Fourier-space bias, b(k), 
defined as the ratio of the power s pectrum for densit y 
peaks and that of dark matter (e.g. Dalai et al.l (2008); 
iMatarrese fc Verdd (|2008l ); IWagner fc Verdd (|2012l )). How- 
ever, an arguably more intuitive measure of the bias is in 
real space, where we can directly obtain information on the 
clustering amplitude of density peaks separat ed by comovin g 
distance r. Following the pioneering work of lKaiserl l|l984) . 
we define 



6 2 (r) = 



(10) 



where r is the comoving length in Eulerian space. The cor- 
relation function f(xi,X2) = (5(xi), <5(xa)) where <5(x) is 
the overdensity field and r — |xi — X2I. Similarly, £ p k is 
the 2-point correlation of density peaks corresponding to 
MHs of mass M. Whilst the real and Fourier space biases 
deal with the same physics of clustering in the overdensity 
field (and they are indeed equivalent in the Gaussian case), 
the real-space bias is related directly to the joint probabil- 
ity distribution of finding two overde nse regions ex ceeding a 
collapse threshold in a given volume (|Kaiserlll984T ). As NG, 
by definition, distorts this probability distribution from the 
Gaussian, the change in the real-space bias seems a natural, 
measurable quantity which can b e calculated, fo r instance, 
using a saddle-point expansion (Valagcas 2010 ) or Edge- 
worth expansion about the Gaussian (IChongchitnan fc Silkl 
2011). The Fourier bias, on the other hand, would be more 
useful when working with the power spectrum and bis pec- 
tra from different NG shape templates iMatsubaraj (I2012T), or 
when redshift-space distortions are incorporated iMao et al.l 
(|2012f) . In this work, however, the real-space bias suffices for 
the calculation of the 21cm fluctuations (|5). 

To date there have only been a handful of calculations 
of the real-space bias in the presence of NG (and unfortu- 
nately b(k) and b(r) are not related via a straightforward 
Fourier tran sform). One such calcu lation is the saddle-point 
approach in IValageasI <|2010l, l2009h . Whilst this formalism 
has been shown to agree with simulations of massive ha- 
los (~ 10 13 AfQ), it has yet to be tested against simula- 
tions of MH-scale resolution. We recognise that there may 
be limitations to the saddle-point formalism in the strongly 
non-linear regime. At the same time, there is not yet any 
convincing analytic model for the non-Gaussian bias in this 
regime, and thus we appeal to the saddle-point approach 
in this work as a first analytical approach to the prob- 
lem. Recent progress in non-Markovian excursion-set theory 
llAdshead et al.ll2012l ; iMusso et al]|2012l ; iParaniape fc Shethl 
120121 ) should soon allow a more accurate calculation of the 
NG bias down to much smaller masses, and high-resolution 
simulations will be needed to elucidate the gas physics on 
such sc ales. 

In (|Chongchitnan fc Silkl 120121 ) . we studied the saddle- 
point approach in detail and found that when b(r) is aver- 



aged over all separation lengths r within the volume that 
the MHs occupy (which in this case is a sphere with radius 
L), the result is the volume- averaged bias, b(M, z), which is 
well-approximated by 



b(M, z) « [1 + (6/5)/ NL A"(z)L 2 ]6 G (M, z), 



(11) 



where 6g is the Gaussian bias and the constant K can be 
determined from calculating b(r) on some fixed scale. K(z) 
roughly grows as a linear function of z as shown in the top 
panel of Fig. ®. 

Using these results and combining it with the flux- 
weighting ([9]), we plot the average bias, j3(z), for /nl = 0, 0.1 
and 1 in the middle panel of Fig. [1] Clearly j3 is sensitive 
to |/nl1 < 1, which boosts the integrand in the nominator 
of (|9} particularly on mass scales around M max - We have 
also checked that the "Gaussian" curves can be closely re- 
produced using the linear bias of iMo fc White! |l996), 111 
agreement with llliev et alJ (|2003T ). 

Our main results are shown in the bottom panel of Fig. 
[T] which shows the redshift variation of the rms 21cm fluc- 
tuations ((8)l for /nl < 1. The peak structure of these curves 
is the result of the competition between terms in Eq. [51 
with a v and 5Tb decreasing with z (see fig. 2 of llliev et al.l 
(2002)), whlist ft increases as shown in the top panel. In 
this figure, we assume an observation bandwidth Av — 1 
MHz and beam angular diameter AO =9 arcminutes. Also 
overlaid are two increas ing curves corresponding to the noise 
(|Furlanetto et alJlaKKSh 
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The noise thresholds assume total effective areas A to t = 10 
m 2 ("LOFAR") and 10 5 m 2 ("SKA"), with integrati o n tim e 
t int = 1000 ho u rs in both cases (see Ide Vos et ail [|20Qg(k 
iDewdnev et al.l l|2009l l for detailed specifications). These 
curves show that even /nl as small as 0.1 will boost the 
fluctuations from a few mKs to tens of mK. Hence, there 
are good prospects of detecting a small NG via the 21cm 
emission from high-redshift minihalos with upcoming radio 
telescopes. 

Moreover, we have checked that the non- Gaussian ef- 
fects on (<5T 2 ) cannot be easily reproduced by changing each 
of the fiducial cosmological parameters (within the observa- 
tional limits). For instance, increasing <rg from 0.8 to 0.9 re- 
sults in < O.lmK increase across the redshift range shown. A 
Fisher matrix analysis could be performed to shed light on 
parameter correlations, but we shall leave this for a future 
investigation. 



4 DISCUSSION 

Here we consider a number of caveats for the results in Fig. 
[1] mainly coming from the fact that MHs are relatively small 
objects whose dynamics are governed by nonlinear physics 
on small scales. 

i) Mass Function: MHs are notoriously difficult to resolve 
in iV-body simulations, requiring at least a ~ 20 Mpc co - 
mov ing box and > 10 1 Jeans-mass par t icles llM eiksin 2011) 



Iliev et all (|2003h : IShapiro et all (|2006l i; ICiardi et all 



(2006) for previ ous attempts). The large simulation in 



llliev et al.l (|2012T l. in particular, was able to resolve down 
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Figure 1. Top: The function K(z) (see Eq. llll l for a range of mass 
scales, M. K(z) has a very weak dependence on M and is inde- 
pendent of /nl- Middle: The volume-averaged bias, /3(z), of mini- 
halos (Eq. [9} with /nl = (solid/black), 0.1 (long dashed/red) 
and 1 (short dashed/blue). Bottom: The fluctuations in the 21cm 
brightness temperature, {5T 2 ) 1 / 2 (Eq. [§), assuming a bandwidth 
Au = 1MHz and angular diameter A8 = 9' for the same range 
of /nl- Overlaid (thin solid/green) are the noise thresholds for 
LOFAR, and SKA-like experiments as discussed in the text. 



Figure 2. Fluctuations in the 21cm brightness temperature, as- 
suming the same telescope specifications as the bottom panel of 
Fig. Q with /nl = 0.1. The "LOFAR" and "SKA" noise curves 
in that figure are also reproduced here (note the linear scales). The 
fiducial model (black/solid) assumes the Press-Schechter mass 
function, MH mass threshold M m i n = Mj and radiation inten- 
sity (at Lyc* frequency) J a = 0. The panels show deviations from 
this model when the mass function (top), M m i n (middle) or J a 
(bottom) is varied. 



to 10 4 Mq MHs and found their abundance t o lie between 
the Pres s-Schechter JPress fc Schechterl 1 1974) and Sheth- 
Tormen IjSheth fc Tormenlll999r ) predictions. We now con- 
sider how the 21cm emission from MHs is affected by the 
choice of mass function. 

In the top panel of Fig. [2] we replot the 21cm fluc- 
tuations in the lower panel of Fig. [T] for /nl = 0.1 using 
the ab ove mass function s along with those of I Tinker et al.l 
(2008) and 1 Warren et all (|2006h (note the linear scale here). 
The Tinker mass function is known to predict n(z) lying 
between t he Press-Schechter and the S heth-Tormen predic- 
tions (e.g. IChongchitnan fc Silkl (|2012l )). The Warren mass 
function was found to be accurate for high-redshift ob- 
.j ects down to 10 7 h ~ x Mr when compared with simulations 
l|Lukic et al.l l2007f l. We see that the Press-Schechter and 
Tinker prescriptions gave similarly high amplitudes of (<5T 6 2 ) 
for 2 < 10, whilst the Warren and Sheth-Tormen prescrip- 
tions prefer lower amplitudes. The trends are reversed for 
2 ~ 20. In any case, the uncertainty in the mass function 
does not appear to affect the detection prospects for LOFAR 
and SKA. 

ii) Uncertainty in M m i n : On MH scales there are a number 
of so-called "gastrophysical" effects which may overwhelm 



the imprints of NG. For instance, iTseliakhovich et al.l (|201ll ) 
showed that for 2 > 10, the velocity of dark matter rel- 
ative to baryons is generally supersonic and thus baryons 
can advect out of dark matter potential, leading to the 
possibility that M m i n > M.j. This may be further exacer- 
bated by feedback mecha nisms such as photoevaporation 
of MHs JTliev et al.l 120051) and shock heating in the IGM 
(|Oh fc Haimanll2003inFurlanetto fc Loebll2004l ). 

The centre panel of Fig. [2] shows the effects on {ST 2 } 
when M m i n increases by a factor of 10, 50 and 100 (with 
/nl = 0.1). The result is a suppression in (5T 2 ) across 
all redshifts (although gastrophysical suppressions are ex- 
pected to be more significant at 2 > 10.) Nevertheless, we see 
that the fluctuation amplitudes are generally robust against 
changes in M m i n . 

Hi) Lya pumping: The 21cm signals from MHs are unlikely 
to be completely immune to the Wouthuysen-Field mech- 
anism, which redistributes the spin states and couples the 
spin temperature to that of radiation sources. This means 
that as the radiation intensity increases, Ts — > Tk and the 
MH emission will be more and mo r e sup pressed. 

Following IChuzhov fc Shapiro! |2006), we introduce the 
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Lya coupling of the form 
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Vo 



1.3 x 10 



_ la f J«T. 1 exp (-0.3(1 + 2 ) 1 / 2 T" 2/3 



A W T K 



1 + 0.47^ 1 



, (12) 



where J a parametrizes the intensity of the radiation sources 
at the Lya frequency (in units of 10 -21 erg cm -2 s -1 Hz -1 
In the bottom panel of Fig. [2] we plot the effects of 



sr 



J a = 0.1,0.5, and 1 on the 21cm fluctuations. Indeed we 
see a strong suppression of (5T 2 ), with J a > 1 capable of 
effectively cancelling the boost from NG. 

However, a more serious issue concerning Lya pump- 
ing is that the MH signals will be completely submerged 
under a huge absorption signal from the IGM , which con- 
tains much more mass tha n that in MHs |Oh fc Mackl 
2003; iFurlanetto fc Obi l200rj ; lYue et al l 120091 1. In our ex- 
ample with J a = 0.1, we find an absorption amplitude of 
|(<5T b 2 )| ~ 20 mK, increasing to 100 mK when J a = 1. We 
therefore conclude that if MHs are indeed exposed to strong 
Lya pumping, their 21cm emission will not be visible unless 
/nl > 1. 



5 CONCLUSIONS 

We have shown that a small amplitude of primordial non- 
Gaussianity may be detectable via the fluctuations in the 
21cm emission of high-redshift MHs. Even with /nl < 1, we 
showed that the strong enhancement in the bias leads to a 
significant increase in the amplitude of the fluctuations, as 
seen in Fig. [T] There are good prospects for such a detection 
by the next generation of large radio telescopes such as the 
SKA. 

Our conclusions rely on a number of assumptions on the 
physics of MHs at high redshift. The analytic formalism used 
to calculate the bias assumes some extrapolations from clus- 
ter scales on which the theory has been well-tested. An im- 
proved calculation awaits a fuller understanding of nonlinear 
collapse in the presence of NG, perhaps with help from the 
resurgence of interests in excursion set theory. The results 
presented here are robust against the assumed mass function 
and minimum MH mass threshold. However, the fluctuations 
are sensitive to the presence of radiation background, since 
the Wouthuysen-Field effect is capable of overwhelming the 
MH signals with the IGM line- absorption signals. 

Nevertheless, our conclusions still apply to MHs that 
are sufficiently isolated from UV sources with no strong 
feedback. This class of MHs provides a new probe for pri- 
mordial non-Gaussianity which has so far been unexplored. 
For MHs that are subjected to strong background radia- 
tion, it may be possible for the MH signals to dominate if 
non-Gaussianity is much larger than the CMB bounds (per- 
haps with /nl ~100). However, this awaits high-resolution 
simulations to elucidate high-redshift gastrophysics, cos- 
mic reionization and the nonlinear bias given highly non- 
Gaussian initial conditions. Needless to say, this will be an 
extremely challenging task. 

We thank Patrick Valageas and, in particular, Ilian Iliev 
for illuminating discussions. We also thank the referee for 
comments that led to major improvements of the first ver- 
sion. SC is grateful for the support of Lincoln College, Ox- 
ford, where part of this work was completed. 
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